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In this paper, we explore the interior dynamics of neutral and charged black holes. Scalar collapses 
in flat, Schwarzschild, and Reissner-Nordstrom geometries are simulated. We examine the dynamics 
in the vicinities of the central singularity of a Schwarzschild black hole and of the inner horizon 
of a Reissner-Nordstrom black hole. In simulating scalar collapses in Schwarzschild and Reissner- 
Nordstrom geometries, Kruskal and Kruskal-like coordinates are used, respectively, with the presence 
of a scalar field being taken into account. It is found that, besides near the inner horizons of Reissner- 
Nordstrom and Kerr black holes, mass inflation also takes place near the central singularity in neutral 
scalar collapse. Approximate analytic expressions for different types of mass inflation are partially 
obtained via a close interplay between numerical and analytical approaches and an examination of 
the connections between Schwarzschild black holes, Reissner-Nordstrom black holes, neutral collapse, 
and charge scattering. We argue that the mass inflations near the central singularity and the inner 
horizon are related to the localness of the dynamics in strong gravity regions. This is in accord with 
the Belinskii, Khalatnikov, and Lifshitz conjecture. 
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I. INTRODUCTION 

The internal structure of black holes and spacetime sin¬ 
gularities have been intriguing and basic topics in gravit¬ 
ation and cosmology, theoretical and realistic [TH^. Ex¬ 
ploring the internal structure of black holes and space- 
time singularities inside black holes can lead to a better 
understanding on black hole physics, gravitation, and 
quantum physics. Because of similarities between the 
singularities inside black holes and the singularity in the 
early Universe, this study may also shed light on cosmo¬ 
logy. 

It is widely believed that in reality, gravitational col¬ 
lapses may produce rotating black holes. Price’s theorem 
states that gravitational radiation, produced on the sur¬ 
face of a collapsing star, carries away all the initial fea¬ 
tures of the star’s gravitational held, except the mass, 
charge, and angular momentum parameters [5]. As a 
next step, there is naturally the question of what the fl- 
nal state of the internal collapses might be. A simpler 
version of the question is how, in reality, inside a rotat¬ 
ing black hole, particles from the accretion disk can affect 
the internal geometry of the black hole. 

On the spacelike singularity curve inside a Schwarz¬ 
schild black hole, the Kretschmann scalar curvature di¬ 
verges. Then the maximal globally hyperbolic region 
deflned by initial data is inextendible. However, inside 
charged (Reissner-Nordstrom) and rotating (Kerr) black 
holes, the central singularity is timelike. The globally 
hyperbolic region is up to the Cauchy horizon, and the 
spacetime is extendible beyond this horizon to a larger 
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manifold. The Reissner-Nordstrom inner (Cauchy) hori¬ 
zon is a surface of infinite blueshift, which in turn may 
cause the inner horizon unstable [7]. This motivated the 
strong cosmic censorship conjecture, which states that 
for generic asymptotically flat initial data, the maximal 
Cauchy development is future inextendible. For math¬ 
ematical explorations of the internal structures of charge 
black holes, see Refs. [iin]. For reviews on the Cauchy 
problem in general relativity and strong cosmic censor¬ 
ship, see Refs. [10] and respectively. 

As the singularities are approached, the tidal force di¬ 
verges, and classical general relativity does not apply. 
It is nevertheless important to explore the dynamics in 
extremely strong-gravity regions in classical general re¬ 
lativity. Belinskii, Khalatnikov, and Lifshitz (BKL) in¬ 
vestigated the asymptotic behaviors in the vicinity of a 
spacelike singularity and found that the generic solution 
has chaotic oscillations of the Kasner axes as the singular¬ 
ity is approached [I2HI4] . When a massless scalar field 
is present, the oscillations will be destroyed, and “the 
collapse is described by monotonic (but anisotropic) con¬ 
traction of space along all directions” [Mil]. Moreover, 
in the vicinity of the singularity, “the variation of the 
gravitational held from one location to the next can be 
neglected — what is much more important is the way 
gravity changes over time” m- The BKL conjecture 
was verified numerically in the singularity formation in a 
closed cosmology in Refs. [Mni], and in the dynamics of 
a test scalar held approaching the singularity of a black 
hole in Ref. [H]. In Ref. [H], the BKL conjecture in the 
Hamiltonian framework was examined under the concern 
of loop quantum gravity. In Ref. [ID], the BKL conjec¬ 
ture was tested in spherical scalar collapse in dark energy 
f{R) gravity in the Einstein frame. Note that in the Ein¬ 
stein frame, in the vicinity of the central singularity, the 
scalar degree of freedom (j)[= (•\/3/21n/')/V87rG] dom- 
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inates the physical scalar field and the potential term for 
(p, where f'=df/dR and R is the Ricci scalar. Therefore, 
the dynamics in this case is essentially the same as that 
in general relativity, and it is expected that the BKL con¬ 
jecture is also valid in spherical scalar collapse in general 
relativity. (In fact, this is also confirmed in the current 
paper.) Compared to such efforts that have been made 
on verifying this conjecture near spacelike singularities, 
little work has been done on examining whether the BKL 
conjecture is valid in the vicinities of the inner horizons 
of charged and rotating black holes. 


A. Mass inflation 

The backreaction of the radiative tail from a grav¬ 
itational collapse on the inner horizon of a Reissner- 
Nordstrom black hole was investigated by Poisson and 
Israel [2T1 [22] . It was shown that due to the divergence 
of the tail’s energy density occurring on the inner hori¬ 
zon, the effective internal gravitational-mass parameter 
becomes unbounded. This phenomenon is usually called 
mass inflation. These arguments were extended to the 
rotating black hole case in Ref. [53] . 

In Refs. [U] [55], approximate analytic expressions 
were obtained by considering a simplified model in which 
the perturbations were modeled by cross-flowing radial 
streams of infalling and outgoing lightlike particles. To 
get more information, some numerical simulations in 
more realistic models have been performed. The dynam¬ 
ics of a spherical, charged black hole perturbed nonlin- 
early by a self-gravitating massless scalar field was nu¬ 
merically studied in Refs. |511129j . Under the influence of 
the scalar field, the inner horizon of a charged black hole 
contracts to zero volume, and the central singularity is 
converted from timelike into spacelike. The mass infla¬ 
tion phenomenon was observed. In Refs. [301131], with 
regular initial data, spherical collapse of a charged scalar 
field was simulated. An apparent horizon was formed. 
A null, weak mass-inflation singularity along the Cauchy 
horizon and a final, spacelike, central singularity were 
obtained. The same process was investigated rigorously 
in Ref. [55]. In addition, gravitational collapses in some 
modified gravity theories have been studied numerically. 
Spherical scalar collapse in f{R) gravity was simulated 
in Ref. [2D]. Asymptotic analysis was implemented in the 
vicinity of the singularity of a formed black hole. Spher¬ 
ical collapse of a neutral scalar field in a given spherical, 
charged black hole in Brans-Dicke theory was investig¬ 
ated in Ref. [33]. Spherical collapses of a charged scalar 
field in dilaton gravity and /(i?) gravity were explored in 
Refs. [51] and [3S], respectively. 

In numerical relativity, it is important to connect ap¬ 
proximate analytic candidate expressions with numerical 
results. In Refs. [361137], the features of the Cauchy ho¬ 
rizon singularity in charge scattering were studied. Ana¬ 
lytic and numerical results were compared at some steps. 

Despite many efforts that have been spent on the mass 


function near the inner horizons of charged and rotating 
black holes, little work has been done in an even sim¬ 
pler case: mass function in the vicinity of the central 
singularity of a Schwarzschild black hole. So far this has 
remained an unexplored area of work. 

In this paper, we use the following notations: 

(i) Neutral collapse: neutral scalar collapse toward a 
black hole formation. 

(ii) Neutral scattering: neutral scalar collapse in a 
(neutral) Schwarzschild geometry. 

(iii) Charge scattering: neutral scalar collapse in a 
(charged) Reissner-Nordstrbm geometry. In this 
process, the scalar field is scattered by the inner 
horizon of a Reissner-Nordstrbm black hole. 


B. New results 

In this paper, we explore neutral collapse, neut¬ 
ral scattering, dynamics in Schwarzschild and Reissner- 
Nordstrbm geometries, and charge scattering. The con¬ 
nections between such processes will also be examined. 

We simulate neutral collapse and investigate the 
asymptotic dynamics in the vicinity of the central sin¬ 
gularity of the formed black hole via mesh refinement. 
Approximate analytic solutions near the central singu¬ 
larity are obtained. We find that, because of the backre¬ 
action of the scalar field on the geometry, mass inflation 
also takes place near the central singularity. In neutral 
scattering, similar results are obtained. 

We explore the dynamics in a Reissner-Nordstrbm geo¬ 
metry and charge scattering. By comparing these two 
processes, we investigate the causes of mass inflation. We 
seek further approximate analytic solutions with the fol¬ 
lowing improvements. Usually, double-null coordinates 
are used in studies of mass inflation in spherical sym¬ 
metry. In the line element of double-null coordinates, 
the two null coordinates u and v are present in the form 
of product dudv. In the equations of motion, mixed de¬ 
rivatives of u and v are present quite often. In this paper, 
we use a slightly modified line element, in which one co¬ 
ordinate is timelike and the rest are spacelike. In this 
case, in the equations of motion, spatial and temporal 
derivatives are usually separated. This simplifies the nu¬ 
merical formalism and helps to obtain approximate ana¬ 
lytic solutions. In addition, we compare numerical res¬ 
ults and approximate analytic solutions closely at each 
step. We compare the dynamics for Schwarzschild black 
holes, Reissner-Nordstrbm black holes, neutral collapse, 
neutral scattering, and charge scattering. With the en¬ 
lightenment of high-resolution numerical results, we treat 
the system as a mathematical dynamical system rather 
than a physical one, examining the contributions from all 
the terms in the equations of motion. 

According to the strength of the scalar field, charge 
scattering can be classified into five types as follows: 
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(i) Type I: spacelike scattering. When the scalar field 
is very strong, the inner horizon can contract to zero 
volume rapidly, and the central singularity becomes 
spacelike. The dynamics near the central singularity 
is similar to that in neutral collapse. 

(ii) Type II: null scattering. When the scalar field is in¬ 
termediate, the inner horizon can contract to a place 
close to the center or reach the center. In this type 
of mass inflation, for each quantity, the spatial and 
temporal derivatives are almost equal. In the case 
of the center being reached, the center is null. This 
type has two stages: early/slow and late/fast. In 
the early stage, the inner horizon contracts slowly, 
and the scalar field also varies slowly. In the late 
stage, the inner horizon contracts quickly, and the 
dynamics is similar to that in the spacelike scatter¬ 
ing case. 

(hi) Type III: critical scattering. This case is on the edge 
between the above two cases. The central singular¬ 
ity becomes null. 

(iv) Type IV: weak scattering. When the scalar field 
is very weak, the inner horizon contracts but not 
much. 

(v) Type V: tiny scalar scattering. When the scalar 
field is very tiny, the influence of the scalar field on 
the internal geometry is negligible. 


In this paper, we will explore the dynamics of the first 
four types of mass inflation, and obtain approximate ana¬ 
lytic solutions for the first two. 

Some similarities between neutral and charged mass in¬ 
flations are obtained. The gravitational/electric field(s) 
is(are) strong in the vicinities of the central spacelike 
singularity of a Schwarzschild black hole and the inner 
horizon of a Reissner-Nordstrom black hole. Therefore, 
as verified by numerical results, the BKL conjecture ap¬ 
plies at both cases. It is found that the BKL conjecture 
can interpret well how the mass inflation happens. The 
dynamics at strong gravity regions is local. Then at such 
regions, the Misner-Sharp mass function does not provide 
global information on the black holes. 

We explore scalar collapses in both general relativity 
and f{R) gravity. The results are similar. For simpli¬ 
city, we focus on scalar collapses in general relativity in 
this paper, and present the results in f{R) gravity in a 
separate paper [55] . 

This paper is organized as follows. In Sec. jnj we build 
the framework for charge scattering, including action for 


charge scattering and the coordinate system. In Sec. Ill 


we set up the numerical formalism for charge scattering. 
In Sec.jTVj neutral collapse will be studied. In Sec.jVj we 
discuss neutral scalar scattering in a Schwarzschild geo¬ 
metry. In Sec. VI the dynamics in a Reissner-Nordstrom 
geometry will be examined. In Sec. lYTTl we explore 
charge scattering, obtaining approximate analytic solu¬ 
tions. In Sec. VIII[ we consider charge scattering with a 


weak scalar field. In Sec. IX the results will be summar¬ 
ized. 

In this paper, we set G = c = 47reo = 1. 


II. FRAMEWORK FOR CHARGE 
SCATTERING 


In this section, we build the framework for charge scat¬ 
tering, in which a self-gravitating massless scalar field col¬ 
lapses in a Reissner-Nordstrom geometry. For compar¬ 
ison and verification considerations, we use Kruskal-like 
coordinates, and set up the initial conditions by modify¬ 
ing those in a Reissner-Nordstrom geometry with a phys¬ 
ical scalar field. 


A. Action 


The action for the charge scattering system in general 
relativity is 

s=j ( 1 ) 


with 




( 2 ) 


Cp 
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(3) 


R/(167rG), and Cp are the Lagrange densities for 
gravity, a physical scalar field "0, and the electric field 
for a Reissner-Nordstrom black hole, respectively. R is 
the Ricci scalar, and G is the Newtonian gravitational 
constant. F^i, is the electromagnetic-field tensor for the 
electric field of a Reissner-Nordstrom black hole. 

The energy-momentum tensor for a massless scalar 
field Ip is 


rp{‘tP) 


2 5{^\C^) 
V\9\ ^9^'' 




The dynamics for the scalar field tp is 
Klein-Gordon equation. 


(4) 

described by the 


Dtp = 0. (5) 

The electric field of a Reissner-Nordstrom black hole 
can be treated as a static electric field of a point charge 
of strength q sittin g at the origin r — 0. In the Reissner- 
Nordstrom metric the only nonvanishing compon¬ 

ents of are Ftr = —Ftr = —qlr^. The corresponding 
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Figure 1: Initial and boundary conditions for charge scat¬ 
tering. Initial slice is at t = 0. Definition domain for x is 
[—Xb Xb]. u = {t — x)/2 and v = {t + a:)/2. 



X 


Figure 2: Contour lines for r defined by Eq. (111 in a Reissner- 
Nordstrom geometry with m = 1 and q — 0.7. Although the 
exact inner horizon is at regions where uv and are 

infinite, r can be very close to the inner horizon r = r- even 
when uv and {t^ — take moderate values. 


energy-momentum tensor for the electric field is |39| 
2 6{y^^\CF) 

= ^ ( 6 ) 
= ^diag(-l,-1,1,1). 

Although Eq. Q is obtained in the Reissner-Nordstrom 
metric, it is valid in any coordinate system, since as seen 
by static observers, the electromagnetic field should be 
purely electric and radial [5^ [32] . We denote the total 
energy-momentum tensor for the source fields as 


^(total) _ I rr\{F) 

^ flU ^ flU ' ^ fiv 


(7) 


B. Coordinate system 

In the studies of mass inflation, the double-null co¬ 
ordinates described by Eq. ^ are usually used, 

ds^ = dudv + r'^dVl^, (8) 


where a and r are functions of the coordinates u and v. u 
and V are outgoing and ingoing characteristics (trajector¬ 
ies of photons), respectively. For convenience, in this pa¬ 
per, we use a slightly modified form described by Eq. , 
obtained by defining u = (t — x)/2 and v = {t + x)/2 [40] . 

ds^ = e~^'^ {—dt^ + dx^) + r^dfl'^. (9) 

This set of coordinates is illustrated in Fig. Similar to 
the Schwarzschild metric, the Reissner-Nordstrom metric 
can be expressed in Kruskal-like coordinates [41] (also see 
Refs. [2Tll22ll39lll2|)- So for ease and intuitiveness, we 
set the initial conditions close to those of the Reissner- 
Nordstrom metric in Kruskal-like coordinates, taking into 
account the presence of a physical scalar field ip. 

In the form of Ref. [32], the Reissner-Nordstrom metric 
in Kruskal-like coordinates in the region of r > r_ can 
be expressed as 

fc_|_ 

ds^ = ^6-2'=+’' - ly^"-' {-de+dx^)+r^dn^ 

( 10 ) 

where r±{= m± m? — q^) and k±[= {r± — r;p]/(2r^)] 
are the locations and surface gravities for the outer and 
inner horizons of a Reissner-Nordstrom black hole, re¬ 
spectively. m and q are the mass and charge of the black 
hole, respectively. r{t,x) is defined implicitly below |43] . 



Some details for this set of coordinates are given in 
Appendix [^ In this set of coordinates, as implied by 
Eq. ( jllj ), the exact inner horizon is at regions where uv 
and (f^ — x'^) are infinite; however, it is found that, even 
when uv and (t^ —x^) take moderate values, r still can be 
very close to the inner horizon, e.g., r = (1 -I- 10“^°)r_. 
(See Fig. i) Therefore, at such places, the interaction 
between the scalar field and the inner horizon still can be 
very strong, and then we can investigate mass inflation 
numerically. 

This formalism has several advantages as follows: 

(i) In the line element Q, one coordinate is timelike 
and the rest are spacelike. This is a conventional 
setup. It is more convenient and more intuitive to 
use this set of coordinates. For the set of coordin¬ 
ates described by Eq. (|^, in the equations of mo¬ 
tion, many terms are mixed derivatives of u and v; 
while for the set of coordinates described by Eq. §, 
in the equations of motion, spatial and temporal de¬ 
rivatives are usually separated. 

(ii) We set initial conditions close to those in the 
Reissner-Nordstrom metric. Consequently, with the 
terms related to the scalar field being removed, we 
can test our code by comparing the numerical res¬ 
ults to the analytic ones in the Reissner-Nordstrom 
case conveniently. Moreover, comparing dynamics 
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for scalar collapse to that in the Reissner-Nordstrom 
case helps us to obtain intuitions on how the scalar 
field affects the geometry. 

(iii) The interactions between a scalar field and the geo¬ 
metry are local effects. In Refs. [211 El], the space 
between the inner and outer horizons are compacti- 
fied into finite space. This overcompactification, at 
least to us, makes it a bit hard to understand the 
dynamics. In the configuration that we choose, the 
space is partially compactified, and the picture of 
charge scattering turns out to be simpler. 


r^yy + 2r,„CT ^ -I- i'Krilj'^y = 0. 


(17) 


Via the definitions oiu= {t — x)l2 and v = {t + x)/2, the 
constraint equations can be expressed in (f, x) coordin¬ 
ates. Equations ( jlTj ) — ( |I^ and ( [l7| ) -I- (16) generate the 
constraint equations for the {tx} and {it} + {xx} com¬ 
ponents, respectively. 


r,tx + r,tcr,x + -|- = 0, 


(18) 


r,tt + r,xx + + r,x<J,x) + 47rr('0ft -h V’^a,) = 0- (19) 


III. NUMERICAL SETUP FOR CHARGE 
SCATTERING 

In this section, we set up the numerical formalisms for 
charge scattering, including field equations, initial con¬ 
ditions, boundary conditions, discretization scheme, and 
tests of numerical codes. 


A. Field equations 

In this subsection, we list the field equations for charge 
scattering. Details for the Einstein tensor and the 
energy-momentum tensor of a scalar field are given in 
Appendix [b| In double-null coordinates ([^, using 


Gl+G% = 8 tt 


2^(total)^ _j_ ^(total)^ 


one obtains the equation of motion for r. 


B. Initial conditions 

We set the initial data to be time symmetric, 

r,t = O’,* = ^/’,t = 0 at t = 0. (20) 

Therefore, in this configuration, the constraint equation 
(181 is satisfied identically. Note that, in this configur¬ 
ation, the values of r,* and cr,* at t = 0 are the same as 
those in the Reissner-Nordstrom metric case. 

We set the initial value for ^ as 


i’\t=o = aexp 


(x - Xo) 


21 


( 21 ) 


where a, b, and xq are parameters. The initial value for 
a is defined to be the same as the corresponding one in 
the Reissner-Nordstrom metric case (10), 


o-2cr| 


lt=0 


= e 


-2a 


iRN 

lt=0 


— - 1 

\r_ 


U- 


fc-l 


/ ( 22 ) 

+ r^xx) ~ ’’ft + = 6 ( 1-^ j , (12) where r is defined by Eq. (111 with < = 0. We obtain 

^ ^ / tViP initial valiip fnr r in plnar 0 ’p Qpattprincr at / = 0 lav 


where r,* = dr/dt and other quantities are defined ana¬ 
logously. Eor simplicity, we define rj = r'^ and integrate 
the equation of motion for rj, instead. The equation of 
motion for rj can be obtained by rewriting Eq. (121 as |40j 


the initial value for r in charge scattering at t = 0 by 
combining Eqs. (12) and (19), 


r,xx = - r^tcr,t - r^x<J,x + 


2r 


- V,tt + V,xx = 2e ( 1 - ^ 


(13) 


- 27rr(t/>2 -hV-f,)-h —e 




(23) 


Gg = 87rr(*°**^9^ provides the equation of motion for tr. 


- cr.tt -h cr,xx + ^ + 47r(V'^j - ^\) + e = 0. 

r ’ ’ G 

( 14 ) 

In double-null coordinates, the equation of motion for 
Ip (1|) becomes 


We set r^x = cr,x = Q at the origin (x = 0,t = 0) as 
in the Reissner-Nordstrom metric case. The definition 
domain for the spatial coordinate x is [—Xb x;,]. Then 
r(x = x,t = 0) can be obtained by integrating Eq. (23) 


- ip,tt + i’.xx + + r,xi’,x) = 0 . 


(15) 


The {mm} and {mm} components of the Einstein equa¬ 
tions yield the constraint equations. 


r,uu + 2r,„CT „ -I- dTrrV'f, = 0, 


(16) 


via the fourth order Runge-Kutta method from x = 0 to 
X = ±Xf,, respectively. 

In this paper, we employ the finite difference method. 
The leapfrog integration scheme is implemented, which 
is a three-level scheme and requires initial data on two 
different time levels. With the initial data at t = 0, we 
compute the data at t = At with a second-order Taylor 
series expansion as done in Ref. |43j . Take pj as an ex¬ 
ample, 

'ip\t=At = V'|t=o + ’/',t|t=oAt -l- -■)/),**|t=o(At)^. (24) 
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Figure 3: Numerical evolution scheme. 


The values of ’4>\t=o and V',t|t=o are set up as discussed 
above, and the value of '4’.tt\t=Q can be obtained from the 
equation of motion for ^ (15). 


Up to this point, the initial conditions are fixed with all 
the field equations being taken into account. The first- 
order time derivatives of r, ct, and ip aA, t = 0 described 
by Eq. (20) ensure that the constraint equation (18) is 


satisfied. The equation for r at t = 0 expressed by 
(23) implies that the constraint equation (19) is satisfied. 


Computations of r, cr, and ip ait = At via a second-order 
Taylor series expansion, as expressed by Eq. (24) for the 
case of Ip, satisfy all the equations of motion. 


C. Boundary conditions 


The values of r, cr, and ip at the boundaries of x = ±Xb 
are obtained via extrapolations. In fact, we are mainly 
concerned with the dynamics around x = 0. Therefore, 
as long as X{, is large enough, the boundary conditions 
will not affect the dynamics that we are interested in. 


D. Discretization scheme 


In this paper, we implement the leapfrog integration 
scheme, which is second-order accurate and nondissip- 
ative. With the demonstration of Fig. and using the 
quantity ip as an example, our discretization scheme is 


expressed below. 


V'.t = 

'^,tt = 

'4’,v = 


'ipup - Ipdii 

2At ’ 

Ipup - 2V'hr + V'dn 

{Aty 

V^ur - '0ul - tpdr + 'fpdl 

4Ax • At 

'4’,t +ip,x- 


'^,xx 


, V'rt - V’lt 

*■- = 

Iplt - 2V'hr + V'rt 
(Ax)2 

= V'.t - V'.x, 


We let the temporal and spatial grid spacings be equal. 
At = Ax. 


E. Tests of numerical code 


It is necessary to test the numerical code before we 
explore the results further. We compare the numerical 
results obtained by the code with the analytic ones for 
the dynamics in a Reissner-Nordstrom geometry, and ex¬ 
amine the convergence of the constraint equations and 
dynamical equations in charge scattering. 

For charge scattering, the dynamics in a Reissner- 
Nordstrom geometry in Kruskal-like coordinates is a spe¬ 
cial case, in which the contributions from the scalar field 
are set to zero. This special case has analytic solutions 
expressed by Eqs. (10) and (II). Therefore, we can test 


our code by comparing the numerical and analytic results 
in the Reissner-Nordstrom metric. Set m = 1, q = 0.7, 
and Ax = At = 10“^. We plot the evolutions of r and 
a along the slice (x = 3xI0“‘^,t = t) in Fig. |^a). As 
shown in Fig.|^a), numerical and analytic results match 
well at an early stage, while at a late stage where gravity 
and electric field become strong, the numerical evolutions 
have a time delay compared to analytic solutions. 

When the numerical results are obtained, we substi¬ 
tute the numerical results into the disretized equations 
of motion and constraint equations, and find that they 
are well satisfied. Moreover, the convergence of the con¬ 
straint equations (18) and (19) is examined. We as¬ 


sume one constraint equation is nth-order convergent: 
residual=0(/i"), where h is the grid size. Therefore, the 
convergence rate of the discretized constraint equations 
can be obtained from the ratio between residuals with 
two different step sizes. 


n = log 2 




O 


(25) 


Our numerical results show that both of the constraint 
equations are about second-order convergent. As a rep¬ 
resentative, we plot the results for the {tx} constraint 
equation (18) in Fig. [^b) for the slice {x = x,t = 0.65). 

Convergence tests via simulations with different grid 
sizes are also implemented [HI 03]. If the numerical 
solution converges, the relation between the numerical 
solution and the real one can be expressed by 


^^real = ^" 


0(/i"), 
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Figure 4: Tests of numerical code for charge scattering, (a) Numerical vs anal3rtic results for a Reissner-Nordstrom black hole. 
m — 1, q = 0.7, and Ax = At = 10“'*. The slice is for {x = 3Ax,t = t). This is a special case of charge scattering with the 
contribution of the scalar field being set to zero. Numerical and analytic results match well at an early stage, while at a later 
stage where the gravity and electric field become stronger, the numerical evolutions have a time delay, compared to analytic 
solutions, (b) Numerical tests for the {tx} constraint equation (181 and the evolution of a on the slice (x = x,t = 0.65). They 
are both second-order convergent. 


where is the numerical solution. Then for step sizes 
equal to h/2 and /i/4, we have 


-^real -F'^ +0 



-^real -F^ +0 



Defining ci = — FS and C 2 = F^ — F^, one obtains 

the convergence rate 


n = logs • (26) 

The convergence tests for rj = r^, a, and 4’ ^re investig¬ 
ated. They are all second-order convergent. As a rep¬ 
resentative, the results for cr are plotted in Fig. Sb) 
for the slice {x = x,t = 0.65). The values of the 
parameters in charge scattering in this section are de¬ 
scribed at the beginning of Sec. [vni We use the spa¬ 
tial range of a; G [—10 10 ] and the grid spacings of 
h = Ax = At = 0 . 02 . 


IV. NEUTRAL SCALAR COLLAPSE 

In Ref. [5D], neutral scalar collapse in f{R) gravity 
was explored. In this section, we consider a simpler case: 
neutral scalar collapse in general relativity. We examine 
the dynamics in the vicinity of the central singularity of 
the formed black hole. Mass inflation in the vicinity of 
the central singularity will be discussed. 


A. Numerical setup 

The numerical setup is a simpler version of the one 
in neutral scalar collapse in f{R) gravity discussed in 
Ref. [20] • The dynamical equations for r, 77 , a, and 4 can 
be obtained by setting the terms related to the electric 
field in the corresponding equations in Sec. [Ill A] to zero, 

r{-r^tt+r,xx)-r‘^t + r% = e~'^'^, (27) 

-r],tt+'n,xx = ‘2e~‘^'^, (28) 


~ -|- a^xx + 


r,tt 


+ 47r(V'^t - 4'^,x) = 0. (29) 


- 4,u + 4,XX + -{-r^t4,t + r^x4,x) = 0 . 


(30) 


In the equation of motion for a ( |29| ), the term {r^tt — 
f',xx)lf can create big errors near the center x = r = Q. 
To avoid such a problem, we use the constraint equation 
(16) alternatively [30]. Defining a new variable g 


g = -2a - ln(-r_u). 


(31) 


one can rewrite Eq. (16) as the equation of motion for g, 
= 47r • — • 4'^^. (32) 


In the numerical integration, once the value of r at the 
advanced level is obtained, the value of a at the current 
level will be computed using Eq. (31). 
































We set the initial data as 


r,tt = r^t = = V',t = 0 at t = 0. 

The initial value of ipir) is defined as 


V'(Olt=o 


a ■ tanh 


(r - ro) 


21 


(33) 


(34) 


with a = 0.1, 6=1, and rg = 5. The local Misner-Sharp 
mass m is defined as [46] 






k.2 


)^i- 


2m 


(35) 


(See Ref. [17] for details on various properties of the 
Misner-Sharp mass/energy in spherical symmetry.) Then 
on the initial slice (x = x, t = 0), the equations for r, m, 
and g are [20ll^ 


r 


,x — 



(36) 


m^r = 47rr^ 




(37) 


g.r = 


(38) 


a black hole, the expansion of the outgoing null geodesics 
orthogonal to the apparent horizon is zero |48j . Then in 
double-null coordinates, on the apparent horizon, there 

is uni 


+ r%)^l -=0. (39) 

Using such a property, the apparent horizon is found and 
plotted in Figs.j^d) and[^e). As shown in Fig. (U^e), the 
radius of the apparent horizon is tah ~ 1-8. Therefore, 
a black hole is formed. 

We plot the Misner-Sharp mass function ( [3^ along the 
slices {x = l,t = t) and (a; = 2.5, t = t) in Fig. [^f), from 
which one can see that, near the central singularity, the 
mass function diverges. 


C. Asymptotic dynamics near the singularity curve 
of a Schwarzschild black hole 

In this section, we are mainly interested in the dynam¬ 
ics near the central singularity of the formed black hole. 
As a preparation, in this subsection, we list some basic 
results on the dynamics near the central singularity of a 
(stationary) Schwarzschild black hole. We consider the 
Schwarzschild metric in Kruskal coordinates, 


Set r = m = g = 0 at the origin (x = 0, t = 0). Then the 
values of r, m, and g on the initial slice {x = x,t = 0) 
can be obtained by integrating Eqs. (36)-(38) from the 
center x = 0 to the outer boundary x = Xb via the fourth- 
order Runge-Kutta method. The values of r, u, and '0 at 
t = At can be obtained with a second-order Taylor series 
expansion, as discussed in Sec. |IIIB[ The value of g at 
t = At can be obtained using Eq. (31). 

The range for the spatial coordinate is defined to be 
X S [0 10]. At the inner boundary x = 0, r is always 
set to zero. The term 2{—r^t4’,t + T,x'4’,x)/f in Eq. (30) 
needs to be regular at x = r = 0. Since r is always 
set to zero at the center, so is r^t- Then we enforce if) 
to satisfy '0,x = 0 at x = 0. The value of g at x = 
0 is obtained via extrapolation. We set up the outer 
boundary conditions via extrapolation. The temporal 
and spatial grid spacings are At — Ax = 0.005. 

The numerical code used here is a simplified version 
of the second-order convergent one for spherical scalar 
collapse in f{R) gravity developed in Ref. [20]. 


B. Black hole formation 


ds^ = 


32m^ 

-e 

r 


{—dt^ + dx^) + 


(40) 


with 


e-2'^ = 


32m^ 


(41) 


As discussed in Appendix [C] in the vicinity of the sin¬ 
gularity curve, the ratios between the spatial and tem¬ 
poral derivatives of r can be expressed by the slope of 
the contour line r = 0, K, 



= \K\, 



(43) 


In the vicinity of the singularity curve, we rewrite t as 
t = tg — where tg is the coordinate time on the singu¬ 
larity curve and ^ <C to- With the spatial coordinate x 
being fixed, a perturbation expansion near the singular¬ 
ity curve directly yields 


When the scalar field is weak enough, the scalar held r « (l6m^tg^) ^ . (44) 

will collapse and then disperse. A hat space is left. How¬ 
ever, when the scalar held carries enough energy, the Moreover, combining Eqs. ( [4l] ) and ( [dd] ), one obtains 
scalar held can collapse into a black hole. In this pa¬ 
per, we are interested in the latter case. In Eig. [^ we 
plot the evolutions of r, cr, and 'ip ^^nd the apparent hori¬ 
zon of the formed black hole. On the apparent horizon of 


-In 

2 


32m3 


ilnrw^ln^. (45) 
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Figure 5: Evolutions in neutral scalar collapse, (a)-(c): evolutions of r, a, and ip. In (a) and (b), the time interval between 
two consecutive slices is lOAt = 0.05. In (c), it is 5At = 0.025. (d) and (e) are for the apparent horizon and the singularity 
curve of the formed black hole, (f) The Misner-Sharp mass function along the slices {x = l,t = t) and {x = 2.5, t = t). 
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As shown in Appendix throughout the whole space- 
time of a Schwarzschild black hole, inside and outside the 
horizon and also near the singularity, there is 

+ r%)^l - (46) 

Namely, for a (stationary) Schwarzschild black hole, the 
Misner-Sharp mass function is equal to the black hole 
mass everywhere. 



Figure 6: Spatial derivative vs temporal derivative near the 
central spacelike singularity in neutral scalar collapse. Point 
A and point B are on one same hypersurface r = const, while 
point C is on another one. At point C, in first-order accuracy, 
r,x ~ {rc — rA)/^x and r,t rs (rs — rc)/At. Since ta = rs 
and the absolute value of the slope of the singularity curve 
dt/dx is less than 1, there is |r, 2 ,/r,t| « |At/Aa:| < 1. AA and 
— |a|f are perpendicular to the lines of r = ri and r = r 2 . a is 
a certain number, such that the magnitude of —|a|r has the 
length as shown in the figure. 


D. Asymptotic dynamics in the vicinity of the 
central singularity of the formed black hole in 
neutral collapse 

In this subsection, we discuss the asymptotic dynam¬ 
ics in the vicinity of the central singularity of the formed 
black hole in neutral collapse, by connecting the reduced 
field equations to numerical results and comparing the 
dynamics of neutral collapse to that in Schwarzschild 
black holes. 

We first explain why near the singularity curve of the 
formed black hole, the ratio between spatial and temporal 
derivatives are defined by the slope of the contour lines 
r = const. We take the quantity r as an example. As 
illustrated in Fig. point A and point B are on one 
same hypersurface r = const, while point C is on another 
one. At point C, in first-order accuracy, ~ {fc — 
ta )/and r^t ~ (xb — xc)/^t- Since the singularity 
curve is spacelike, the absolute value of the slope of the 


singularity curve K=dt/dx is less than 1. Consequently, 
there is 



At 

Aa; 


|A|<1. 


Similar arguments yield 


(47) 


,xx 


At 

x,tt 


Ax 




(48) 


This can also be interpreted in the following way. As 
shown in Fig. in double-null coordinates, the time vec¬ 
tor i is not normal to the hypersurface of r = const. 
Then the derivatives in the radial direction have nonzero 
projections on both hypersurfaces of x = const and 
t = const. Along an arbitrary slice {x = const, t = t), 
near the singularity, the ratios between first- and second- 
order spatial and temporal derivatives are equal to |Ar| 
and K^, respectively. 

The quantities a and tp have similar features as de¬ 
scribed by Eqs. (|47| and (481. This can be explained as 


follows. Take the scalar field ip as an example. With the 
illustration of Fig. as this scalar field moves toward 
the center, two neighboring points on this scalar wave 
ip should take close values when they respectively cross 
points C and D (which are on one same hypersurface 
r = const) at two consecutive moments, because these 
two points on the scalar wave are neighbors and the “dis¬ 
tances” AD and BC are more important for their values 
than the difference between these two neighboring points. 
In other words, in the vicinity of the singularity curve, 
gravity is more important than the difference between 
neighboring points on the scalar wave. Alternatively, to 
a large extent, later evolutions in a strong gravitational 
field largely erase away the initial information on the con¬ 
nections between neighboring points. These arguments 
are supported by numerical results. Near the singular¬ 
ity, the evolution of ip is described by ip « Cln^, where 
^ = to — t<Cto and to is the coordinate time on the sin¬ 
gularity curve. In Fig. ^ means AD and BC. The 
parameter C changes slowly along the singularity curve, 
compared to the dramatic running of In f |2(lj . In other 
words, ipA ~ ipB and ipc ~ ipD- Then as in the case of r, 
there is 


i’,x 


At 

'Pd 


Ax 


(49) 


Taking the slice (x = 2, t = t) as an example, we ex¬ 
amine the dynamics along this slice via mesh refinement 
that was implemented in Refs. uniiin] and plot the res¬ 
ults in Fig. As shown in Fig. |^f), in the vicinity of 
the singularity along this slice, the ratios between spatial 
and temporal derivatives are 


X.x ^ XT x ^ '4^,x ^ I X^xx ^ I.XX 

x,t cr^t 'ip.t V V '^4* 



( 50 ) 
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As a comparison, the slope of the singularity curve, which 
is plotted in Fig. [^d), at a; = 2 is if « 0.25. 

Because of the fact that in the vicinity of the singular¬ 
ity, the ratios between spatial and temporal derivatives 
for r, CT, and '0 take similar values, and with the numerical 
results plotted in Figs. M and [T^ the field equations 
for r (27), a (29), and 0 (30) can be simplified as follows: 



rr^tt « -r% 

(51) 


X tt . ,2 

cr.tt ~ -l-47r0(, 

r 

(52) 


2 

r 

(53) 

The asymptotic solutions to Eqs. (51)-(53) are [20] 


r 


(54) 

a - 

s Rln^ -1- (To ~ [;d(l - P) - 47r(7^] In^ -|- (Tq, 

(55) 

0 s 

s (Fine 

(56) 


Substituting Eq. ( [54| ) into Eq. ( [^ yields 


Then we have 



(57) 


E. Causes of mass inflation in neutral scalar 
collapse 


In neutral scalar collapse, the equation of motion for 
a is 


~ ^,xx + 


r,u - r^xx 


+ 47r(02( - 0^,) = 0. (58) 


Then as discussed in the above subsection, the term 
47r(0^j — ip\) is positive and can make a{x,t) greater 
than the corresponding value in the Schwarzschild black 
hole case. [See Eq. (55).] This makes the mass function 
divergent near the singularity as will be discussed below. 

Near the singularity, using Eqs. (54|-(57), there is 




(59) 


Then with Eq. (35), the mass function can be written as 


m= -[1 + e “"(rt-r,,)] 


-(1 - 




(60) 


_ ^2^^3+167rC^g2(To 


-leTrC-" 


In the Schwarzschild black hole case, (7 = 0. As shown 
in Appendix]^ the mass function is always constant and 
equal to the black hole mass. In neutral collapse, the 
parameter /3 does not change much and remains about 
1/2. However, the parameter C is not zero. Then the 
metric quantity a is modified. As a result, the delicate 
balance between r and + r\) is broken. Con¬ 

sequently, as implied in Eq. ( p0| ), near the singularity, the 
mass function diverges: mass inflation occurs. With nu¬ 
merical simulations, the divergence of m near the central 
singularity was also reported in Ref. [35], while further 
explorations were absent. 


As shown in Fig. [5];b), at large- and small-a: regions, 
a approaches negative and positive infinities, respect¬ 
ively. We interpret this difference using Eqs. (55)-(57l, 
Fig.j^c), and three sample sets of results at x = 1, a; = 2, 
and X = 2.5, which are obtained via mesh refinement and 
are plotted in Figs, [^and EE and [To] respectively. As 
shown in Fig. [^c), at large-x regions, the scalar field 
is very weak. This is also illustrated in Figs. I^c) and 
[Toj[b), which show the contributions of all the terms in 
the equation of motion for a at the slices {x = 2,t = t) 
and (x = 2.5, t = t), respectively. From these two figures, 
one can see that at large-x regions, in the vicinity of the 
singularity, in the equation of motion for cr, the contri¬ 
bution of the scalar field 0 is much less than those from 
metric quantities, 47r0^( ^ ~ As shown 

in Figs. I^d) and [Toj[c), on the slices {x = 2,t = t) and 
{x = 2.5, t = t),C = -0.0606 and C = -0.02786. Then 
<1/4. Consequently, as implied in Eq. (55) and as 
plotted in Figs, [^d) and 10 c), in such a case, a is neg¬ 
ative. We also note that because of the negativeness of 
a, the contribution of the term 26“^°^ in the equation of 
motion for r] ( |^ is important. [See Figs.j^b) and[To|[a).] 
This is quite different from the positive case, as plotted 
in Fig. gb). Because of weakness of the scalar field at 
large-x regions, a is not modified very much, compared to 
the corresponding value in the Schwarzschild black hole 
case. As a result, as implied in Eq. ( |60| , the mass func¬ 
tion does not grow as fast as in the strong scalar field 
[See Figs.[§;d),|^e), and[^d).] 


case. 


At small-x regions, before the center x = 0 converts 
from timelike into spacelike, the scalar field is reflected 
at X = 0. In addition, new pulses of the scalar field 
from large-x regions keep coming to the central region. 
Then some of the energy of the scalar field is accumu¬ 
lated around the center. Then the scalar field becomes 
very strong in this region. [See Fig. [^c).] In this case, 
the contribution from the scalar field can be comparable 
to or even greater than those from the metric quantit¬ 
ies, 47r0^( « W tt\ ~ Vtt/r]. (See Fig. m) As shown in 
Fig-i c), C ~ —0.18 in the vicinity oT the singularity 
on the slice (x = l,t = t). There is 47r(7^ « 0.41 > 1/4. 
Consequently, tr grows to be positive. [See Fig.|^b).] Be¬ 
cause of the positiveness of ct, the contributi on from the 
term 2e~^’^ in the equation of motion for rj ( |2^ is small 
compared to those from rj^tt and r] xx- [See Fig. mb).] Be¬ 
cause of strongness of the scalar field, a is significantly 
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Figure 7: (color online). Dynamics in the vicinity of the singularity in neutral collapse on the slice {x = l,t — t). The numerical 
results show that in the vicinity of the singularity, the field equations for r (271, r} (281, a (291, and (30l can be reduced to 
the following formats, respectively: (a) rr^tt ~ 


(b) V.i 


V,x 


(c) cr,tt — r^tt/r « 47rV>^and (d) ~ —2r,t'!/’,t/r. The 


big numerical errors near r = 0 in (a) come from converting Eq. (131 into Eq. (121. 


modified. As a result. 


ically. [See Eq. (60) and Fig. [^d) 


F. Differential form of the mass function near the 
central singularity 

As a further step of exploring mass inflation, we ex¬ 
amine the differential form of the mass function. Using 


Eq. (46) and Einstein equations, one obtains nans], 


£ = ,61) 

where x°' is a coordinate, t or x. {=Tl+T^) is zero for 


the r nassle ss scalar field with V(ip) = 0. [See Eqs. (B9) 
and (BIO).] We denote f, t, and x as unit vectors along 


the radial, temporal, and spatial directions, respectively. 
As illustrated in Fig.[^ f is perpendicul ar to the c ontou r 
lines r = const. Combining Eqs. (61), (B9), and (BIO), 
one obtains the gradient of m, 

dm - dm ^ 
yrn = —t -h -^x 
dt dx 


= -dTrr^ • + V'^,„)[r tf -h r 

= -47rr^ • + K^(-f) 

AA ^ ’ 


(62) 


As implied in Fig. 


AA = 


At ■ Ax 


At 


y/(Af)2 -h (Ax)2 y'l + K^' 


(63) 
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(a) 





Figure 8: Solutions in the vicinity of the central singularity in neutral scalar collapse on the slice (x = l,t = t). (a) Inr = 
aln( + b, a = 0.49948 ± 0.00004, b = 0.2172 ± 0.0005. (b) cr = aln^ + 6, a = -0.1651 ± 0.0001, b = -0.083 ± 0.001. (c) 
= aln^ + b, a = -0.18129 ± 0.00003, b = 0.1322 ± 0.0003. (d) Inm = alnr + b, a = -1.670 ± 0.002, b = -1.258 ± 0.009. 


Then using 


G. Discussions 


Ar Ar At 
AA ^ Af ' Aa' 




(64) 


we have 


^ = ( 65 ) 


AA Ar 


Equation ( |65[ ) implies that in the vicinity of the central 
singularity, the effective density for the scalar field ip is 


Peft = (66) 


which is consistent with the integration form of the mass 
function (60). 


1. In which cases is the mass function equal to the 
black hole mass? 

For stationary black holes, the Misner-Sharp mass 
function is always equal to the black hole mass. This 
is verified for the Schwarzschild black hole case in Ap- 
pendix [C]a nd for the Reissner-Nordstrom black hole case 
in Sec. Ivil 

In spherical symmetry, at spatial infinity, the mass 
function describes the total energy/mass of an asymp¬ 
totically flat spacetime m- In a gravitational collapse 
case, it means the total mass of the collapsing system. 

2. In which cases does the mass funetion not describe 
the black hole mass? 

In the vicinities of the central singularity of a Schwar¬ 
zschild black hole and the inner horizon of a Reissner- 
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Figure 9: (color online). Dynamics and solutions in the vicinity of the central singularity in neutral scalar collapse on the slice 
{x = 2,t = t). (a)-(c): dynamical equations for r, r], and a. (d) a = aln^ + b, a — 0.1980 ± 0.0004, b = —0.112 ± 0.003. 
^ = alnC + fo, a = -0.0606 ± 0.0001, b = 0.2899 ± 0.0008. (e) Inm = alnr + fe, a = -0.1993 ± 0.0006, b = 0.013 ± 0.002. (f) 
Ratios between spatial and temporal derivatives for r, a, and ip. 
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(c) 



10 ’^ 10 ’'' 


r 

(d) 


Figure 10: (color online). Dynamics and solutions in the vicinity of the central singularity in neutral scalar collapse on the slice 
{x = 2.5, t = t). (a) and (b): dynamical equations for rj and a. (c) a = aln^ + b, a — 0.2476 ± 0.0002, b = —0.337 ± 0.002. 
i> = aln^ + b, a =-0.02786 ±0.00003, 6 = 0.2732 ± 0.0003. (d) lnm = alnr±6, a =-0.046 ± 0.001, 6 = 0.072 ± 0.003. 


Nordstrom or Kerr black hole, the dynamics and 
some quantities are local. The mass function is just a 
parameter which varies at each point, not giving global 
information on the black hole mass. 

3. How does one interpret the behaviors of the mass 
function at different circumstances? 

One may interpret this issue via an analog to Newto¬ 
nian gravity. In Newtonian gravity, suppose we want to 
measure the mass M of a source sphere with a point¬ 
like test mass m. Denote the distance between the two 
masses by r. Without perturbations, the gravitational 
force between the two masses is simply F = GMm/r'^. 
However, if another mass 6M passes by the test mass, 
the gravitational field near the test mass can become very 
strong and local, and it is not able to provide accurate 
information on the source sphere. When the perturba¬ 
tion mass is gone or sticks to the source sphere, one can 


measure the mass of the source object accurately again. 

Similarly, in gravitational collapse in general relativity, 
in the vicinity of the singularity of the formed black hole, 
because of the strong interaction between the scalar field 
and the geometry, the dynamics becomes local. Then the 
computed mass function does not provide global inform¬ 
ation on the mass of the black hole or of the collapsing 
system. As discussed in Sec. |IVE[ at large-a; regions, 
most of the energy of the scalar field has been absorbed 
into the black hole, and the tail of the scalar field is very 
weak. As a result, although in this case the mass func¬ 
tion does not seem to approach a fixed value, it grows 
very slowly. [See Fig. 10 [d).] 

4 . Is there a suitable mass definition inside black 
holes? 


Although quite several local definitions of mass have 
been constructed for regular spacetime, little work has 
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been done to provide more alternatives for global con¬ 
servation laws. Considering the large amount of efforts 
that have been spent on formulating proper quasi-local 
mass, it can be even more challenging to develop the “ul¬ 
timate” mass definition inside black holes. The role of 
unavoidable spacetime singularities has also an effect in 
all of these definitions. 


V. NEUTRAL SCALAR SCATTERING 


In this section, we consider neutral scattering, in which 
a neutral scalar collapses in a Schwarzschild geometry. 
The numerical formalism is a simpler version of the one 
in charge scattering that has been constructed in Sec. m 
and it can be obtained by removing the electric terms in 
the held equations presented in Sec. Ill A and replacing 
the Reissner-Nordstrom geometry with a Schwarzschild 
one. 

The parameters are set as follows: 


VI. DYNAMICS IN THE 
REISSNER-NORDSTROM GEOMETRY 


As one more preparation for charge scattering, we dis¬ 
cuss the dynamics in the Reissner-Nordstrom geometry 
in the patch r > r_ in Kruskal-like coordinates 

fc_|_ 

ds^ = {-dt^+dx^)+r^dn^ 

kj^r^ \r_ J 

(67) 

where r and a are defined by 



2cr 2kj. 


e = 




— - 1 
r_ 


i+- 


fc-i 


(69) 


The equations of motion for r and cr are 


(i) Schwarzschild geometry: m = 1. 

(ii) Physical scalar held: 

= aexp [-(x - xo)'^/b], 
a = 0.2, b = 1, and Xq = I. 


'’(-'T.tt + -b 


rl = 



(70) 


~ CT tt + a XX + 


,XX 


= 0 . 

r4 


(71) 


(hi) Grid. Spatial range: x € [—10 10]. Grid spacings: 
Ax = At = 0.005. 


The numerical results for neutral scattering are plot¬ 
ted in Fig. 11 As shown in Fig. ©b), a is always 
negative. Our interpretations for this are the follow¬ 
ing. In neutral collapse, before the center x = 0 changes 
from timelike into spacelike, because of the rehection at 
X = 0 and the consecutive incoming of pulses of the 
scalar held from large x, some of the energy of the scalar 
held is accumulated at the central region. Then the 
scalar held becomes very strong and its contribution can 
be comparable to those from the metric quantities, i.e., 
^ \r^tt/r\ ^ As a result, a can be positive. 

However, in neutral scattering, there is a given Schwarz¬ 
schild geometry, and the center is a spacelike singularity 
from the beginning. So the energy of the scalar held 
simply moves to the central singularity rather than be¬ 
ing rehected. Then the energy of the scalar held has no 
opportunity of accumulating at the central region. Con¬ 
sequently, the scalar held is always weak compared to the 
metric quantities, and cr remains negative. 

We explore the dynamics in the vicinity of the central 
singularity, including the behaviors of the held equations, 
the metric quantities, the scalar held, and the mass func¬ 
tion. The results are similar to those in the weak scalar 
held case in neutral collapse as plotted in Figs. and 
[T0| As a representative, we plot the Misner-Sharp mass 
function along the slice (x = 2, t = t) in Fig. ©f)- 


The mass and charge of the Reissner-Nordstrom black 
hole are set as to = 1 and q = 0.7, respectively. We 
take the slice (x = 0.5, t = t) as a sample slice. The 
dynamics for r and a and the evolutions of r, a, and 
TO on this slice are plotted in Fig. 12 As shown in 
Fig. 12 in Eqs. (70) and ( [7l| ), around the outer hori¬ 
zon, the charge terms are negligible compared to other 
terms. The dynamics is close to that in the Schwarzschild 
black hole case, a is negative and jcrj keeps decreasing 
until r is close to the inner horizon. Then jcrj deceler¬ 
ates and turns around, a approaches positive infinity, 
and approach zero, and r asymptotes to r_. The 
term e ^‘^(—+ r\) asymptotes to zero, as plotted of in 
Fig. [^f). Note that the Misner-Sharp mass function a 
charged black hole is defined as 




2m 


.,,,.=e-(-r:i+r:L)^l- —+ ^. 


(72) 


Then as shown in Fig. [I^f), the mass function remains 
equal to the mass of the Reissner-Nordstrom black hole 
as the inner horizon is approached as expected. 

In the next section, we will show that in charge scatter¬ 
ing, near the inner horizon, due to the additional contri¬ 
bution from a scalar field, a approaches positive infinity 
faster than in the Reissner-Nordstrom black hole case. 
Gonsequently, r decreases faster and then is able to cross 
the inner horizon. As a result, the term + r\) 

and the mass parameter diverge: mass inflation takes 
place. 
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X 

(e) 



Figure 11: Evolutions in neutral scattering, (a)-(c): evolutions of r, a, and The time interval between two consecutive 
slices is 20At = 0.1. (d) and (e) are for the apparent horizon and the singularity curve of the formed black hole, (f) The 
Misner-Sharp mass function along the slice {x = 2, t = t). 
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(d) 



(f) 


Figure 12: (color online). Dynamics in a Reissner-Nordstrom geometry on the slice (x = 0.5, t = t). (a) and (b): dynamical 
equation for r. (c) and (d): dynamical equation for a. (e) and (f): evolutions of r, a, and m. At the early stage, r and a 
move toward zero and negative infinity, respectively, as in the Schwarzschild metric. Later on, due to the repulsive force from 
2e“^'^g^/r^, r decelerates and asymptotes to r_. [See Figs, (a) and (e).] Mainly because of jr'^, a switches the sign and 

approaches positive infinity. [See Figs, (c) and (e).] As the inner horizon is approached, the term asymptotes 

to zero, and the mass function m remains equal to the mass of the Reissner-Nordstrom black hole. [See Fig. (f).] 
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X X 

(e) (f) 

Figure 13: Evolutions in charge scattering. In (a), (b), and (d), the time interval between two neighboring slices is 30At = 0.15. 
In (c), it is 15At = 0.075. (e) and (f) are for the horizons. When the scalar field is strong enough (around x = 2), the inner 
horizon can be pushed to the center, and the central singularity becomes spacelike. When the scalar held is weak enough (e.g., 
—4 < X < —1), the inner horizon does not change much, and the central singularity remains timelike. At the intermediate 
state (e.g., 0 < x < 1.5), the inner horizon contracts to zero, and the central singularity becomes null. The results for the inner 
horizon, especially for x > 2, are not that accurate. We are aware that the inner horizon is actually at inhnity, while r still can 
be very close to r_ when x and t take moderate values. 
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VII. RESULTS FOR CHARGE SCATTERING 

In this section, we explore charge scattering: scalar 
field collapse in a Reissner-Nordstrom geometry. We 
study the evolutions of the metric components and scalar 
field, and obtain approximate analytic solutions. We 
closely compare the dynamics in Schwarzschild black 
holes, Reissner-Nordstrom black holes, neutral scalar col¬ 
lapse, and charge scattering. In addition, in obtaining 
the analytic solutions, we compare the analytic candid¬ 
ate solutions to the numerical results at each step. 

In this section, the parameters are set as follows: 

(i) Reissner-Nordstrom geometry: m = 1, and q = 0.7. 

(ii) Physical scalar field: 

= aexp [-{x - x^Y/b], 
a = 0.08, 6=1, and Xq = 4. 


(iii) Grid. In Sec. VIIA we use the spatial range of 
X G [—10 10] and grid spacings of Ax = At = 0.005. 
In Secs. |VII B|VilE[ we use the spatial range of 
X G [—5 5] and grid spacings of Ax = At = 0.0025. 


A. Evolutions 


(ii) Type II: null scattering. When the scalar field is 
intermediate, the inner horizon can contract to a 
place close to the center or reach the center. In this 
type of mass inflation, for each quantity, the spatial 
and temporal derivatives are almost equal. In the 
case of the center being reached, the central singu¬ 
larity becomes null. This type of scattering has two 
stages: early/slow and late/fast. In the early stage, 
the inner horizon contracts slowly; while in the late 
stage, the inner horizon contracts rapidly. Sample 
slice: {x = 0.5, t = t) in Fig. [U 
and lylml 

(iii) Type III: critical scattering. This case is on the 
edge between the above two cases. When the central 
singularity is reached, it becomes null. Sample slice: 


See Secs. VIIC 


{x = 1.53, t = t) in Fig. 13 See Sec. VIIE 


(iv) Type IV: weak scattering. When the scalar field 
is very weak, the inner horizon does not contract 
much. Sample case: Figs. and See Sec. [VIII| 

(v) Type V: tiny scalar scattering. When the scalar 
field is very tiny, the influence of the scalar field 
on the geometry is negligible. Sample slice: {x = 
—3, t = t) in Fig. 13 


In this paper, we will discuss the first three types. 


1. Outline 


2. Causes of mass inflation and evolutions of r, a, and ‘ip 


In this subsection, we describe the evolutions of r, a, 
and tp that are plotted in Fig. |13[ Examining the equa¬ 
tions of motion (12)-(15), one can see that in the charge 
scattering dynamical system, there are three types of 
quantities as follows: 

(i) Metric components: r and a. They contribute as 
gravity. 


In a Reissner-Nordstrom geometry, in Kruskal-like co¬ 
ordinates expressed by Eq. (67), near the inner horizon. 


although (T asymptotes to positive infinity, {r\ — r\) is 


Then as implied in Eq. (72), m takes a 


much less than Consequently, ap¬ 

proaches zero, 
finite value. 

In charge scattering, the equations of motion for 77=r^ 
and a are 


(ii) Scalar field: ip. It contributes as a self-gravitating 
field. 


V,tt + V,xx = 2e ( 1 - ^ 


(73) 


(iii) Electric field. 
Eq. pi). 


It acts as a repulsive force. See 


Eurthermore, these quantities can be separated into two 
sides: the gravitating side (r, a, and ip) and the repulsive 
side (electric field). The dynamics in charge scattering 
consists mainly of how these quantities interact, includ¬ 
ing how the gravitating and repulsive sides compete. 

According to the strength of the scalar field, charge 
scattering can be classified as follows: 

(i) Type I: spacelike scattering. When the scalar field 
is very strong, the inner horizon can contract to zero 
volume rapidly, and the central singularity becomes 
spacelike. Sample slice: {x = 1.67, t = t) in Eig. 13 
See Sec. Nim 


— CT.tt + O' XX + 


,xx 


■ 47r(V'^4 - ipl 




= 0 . 

(74) 
is 
as 


At the beginning, under our initial condition s, \ip^ 
less than However, as discussed in Sec. 

r decreases toward the central singularity, gravity be¬ 
comes stronger, and \ip^t\ becomes greater than \ip^x\- 
(Also see Figs. ppTl ) Then in this case, the repuls¬ 
ive “force, ”47r(^pp*j,) -|-is greater than the 
corresponding one, in a Reissner-Nordstrom 

geometry. This makes a accelerate faster than in the 
Reissner-Nordstrom geometry. Consequently, the re¬ 
pulsive force from 2e“^'^(q^/r^ — 1) for r]=r'^ is much 
weaker than the corresponding value in the Reissner- 
Nordstrom geometry. As a result, near the inner horizon. 
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\r^t\ is much greater than the corresponding one in the 
Reissner-Nordstrom case. Then {r\ — r\) moves from 
extremely tiny values in the Reissner-Nordstrom metric 
case to moderate values, and r crosses the inner horizon 
r = r- for the given Reissner-Nordstrom geometry. [See 
Figs. 12 [e) and [l5j[a) .] This makes the mass parameter 


diverge: mass inflation occurs. [See Eq. (72).] In other 
words, regarding the causes of mass inflation in charge 
scattering, the scalar held’s backreaction on r is more 
important than that on a. The evolutions of r, cr, and if) 
are plotted in Figs. [T^a )@d). 


negligible, compared to other terms in the equa¬ 
tions. As a result, in the vicinity of the central 
singularity, the dynamics is similar to that in neut¬ 
ral, strong scalar collapse as expressed by Eqs. (51 )- 
(53). [See Eigs. [^a)-[^d) .] Then the quantities 
r, cr, ip, and m take similar forms as those in neut¬ 
ral collapse and neutral scattering as expressed by 

and[T^f).] 


Eqs. (54)-(57) and (60). [See Figs. 16 


C. The late/fast stage of null scattering 


3. Locations of horizons 


In charge scattering, on the horizon, there is 




2m 


> = + rfj^l-+ % = 0. (75) 


We locate the outer and inner horizons using this prop¬ 
erty. The results are plotted in Figs. 13 [e) and [l^f). 


Because of the absorption of the physical scalar field ip, 
the outer horizon increases from the original value of 1.7 
to 2.9. Note that the results for the inner horizon, espe¬ 
cially at regions where x > 2, are not that accurate. We 
are aware that the inner horizon is actually at inhnity, 
while r still can be very close to r_ even when x and t 
take moderate values. 


B. Spacelike scattering 

In this configuration, the scalar field is so strong, such 
that the inner horizon can contract to zero volume rap¬ 
idly, and the central singularity converts from timelike 
into spacelike. Taking the slice {x = 1.67, t = t) as a 
sample slice, we plot the terms in the field equations for 
r, (T, and ip in Fig. and also the evolutions of r, a. 


Ip, and TO in Fig. 15 We investigate the dynamics in the 
vicinity of the central singularity via mesh refinement, 
and plot the results in Fig. 

The strongness of the scalar held causes several con¬ 
sequences as below. 

(i) The quantity r does not decelerate much when it 
crosses the inner horizon of the given Reissner- 
Nordstrom black hole, and it can approach the cen¬ 
ter. At the same time, the mass function grows 


dramatically. [See Fig. 15 


(ii) The central singularity converts from timelike into 
spacelike. 


(iii) The dynamics in spacelike scattering is similar to 
that in neutral, strong scalar collapse. The quant¬ 
ity a takes large positive values, such that in the vi¬ 
cinity of the central singularity, the term e~‘^^cp' jr'^ 
in the equation of motion for r ( 12 ) and t he t erm 
g- 2 (T^ 2^^4 equation of motion for cr (14) are 


When the scalar held is less strong, the inner horizon 
may still contract to zero. However, in this case, the 
central singularity becomes null rather than spacelike. 
The equations of motion remain null: they have similar 
forms as free wave equations, e.g., V'.tt ~ ip,xx- 

In the Reissner-Nordstrom black hole case, near the 
center, the repulsive (electric) force dominates gravity, 
and the central singularity is timelike. In spacelike scat¬ 
tering as discussed in the last subsection, gravity from 
the scalar held and the background geometry dominates 
the repulsive force. As a result, the central singularity 
is spacelike. At the late stage of the null scattering that 
will be studied in this subsection, the scalar held is less 
strong, and the central singularity is null. Because of 
this, one may say that the null scattering is a critical 
case of the competition between repulsive and gravita¬ 
tional forces, in which case the two types of forces have 
a balance. 

Take the slice {x = 0.5, t = t) as a sample slice, we 
plot the terms in the held equations for r, a, and ip in 
and the evolutions of r, a, ip, to, and |1 — in 
, We investigate the dynamics in the vicinity of 


17 


18 


Fig. 

Fig. 

the central singularity via mesh rehnement and plot the 
results in Fig. 

The null scattering has two stages: early/slow and 
late/fast. As shown in Figs. 17 and right after the 
strong collision between the scalar held and the inner ho¬ 
rizon, because of the repulsive force from the electric held 
and the tension force from spatial derivatives, r, a, and 
Ip change slowly. As a result, the mass function to also 
grows slowly. We call this stage the early/slow stage. 
Later on, as r approaches zero, maybe because of strong 
gravity from the central singularity as in the spacelike 
scattering case, these quantities change faster. We call 
this stage the late/fast stage. 

As shown in Fig. Ill when r is very small, the equations 
of motion for r (12), 77 (13), a (14), and ip (15) can be 
rewritten as 


- rrtt 


-rr 

(76) 

~ t 

(77) 

i 47r7):2 ~ 47r7/>2^, 

(78) 























Terms in E.o.M. for ib Terms in E.o.M. for cr Terms in E.o.M. for r 
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Figure 14: (color online). Dynamics along the slice {x = 1.67, t = t) in spacelike charge scattering, (a) and (b), (c) and (d), (e) 
and (f): dynamical equations for r, a, and tp. In this conhguration, the scalar field is so strong, such that r does not decelerate 
much when it crosses the inner horizon of the given Reissner-Nordstrom geometry. 
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Figure 15: (color online). Evolutions along the slice (x = 1.67, t = t) in spacelike charge scattering, (a) evolutions of r and a. 
(b) evolutions of ip and m. The mass function remains equal to the mass of the original Reissner-Nordstrom black hole, m — 1, 
until t « 1.9, and by then mass inflation takes place, (c) ratios between spatial and temporal derivatives for r, a, and ip. When 
t > 1.9, there is \ip,x/ip,t\ ~ \r,x/r,t\. 


2 2 

i’,tt « 'ip, XX « --J'.tV'.t « --r,x'ip,x- (79) 


D. The early/slow stage of null scattering 


Since the above four equations have some similarities to 
the corresponding ones in spacelike scattering, one may 
guess that the quantities r, tr, ip, and m may have ex¬ 
pressions similar to those in spacelike scattering. In fact, 
this guess is verified by the numerical results plotted in 
Fig.lH 

In the spacelike scattering case, near the central sin¬ 
gularity, the ratio \r^xlf,t\ approaches a fixed value—the 
slope of the singularity curve, which is different from 1. 
Therefore, we can obtain the approximate analytic ex¬ 
pression for {r\ — r\) and then for the mass function m. 
However, in the null scattering case, |1 — K'^ \ asymptotes 
to zero. We plot the numerical results of |1 — in 
Fig.[^b), while we do not derive the approximate ana¬ 
lytic expression for |1 — K‘^\. The good thing is that in 
the diverging mass function, compared to the factor e^'^, 
the term |1 — K‘^\ is a minor one. 


As shown in Fig. 17 at the early/slow stage of null 


scattering, the equations of motion for r (12), cr (14), 
and '0 (15) can be rewritten as follows: 


2^2 

(80) 


(81) 

,xx: 

(82) 


The above equations are like free scalar wave equations 
in flat spacetime. The derivatives of one variable (r, tr, 
and Ip) are independent from the derivatives of another. 
Arbitrary functions of {t -I- x) or {t — x) can satisfy the 
above equations, and in principle, the initial conditions 
right after the collision between the scalar field and the 
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(e) (f) 


Figure 16: (color online). Dynamics and solutions in the vicinity of the spacelike central singularity on the slice {x = 1.67, t = t) 
in spacelike charge scattering. In the equations of motion, in the vicinity of the singularity, the terms related to the electric 
field are negligible. Therefore, the equations of motion behave similarly to those in neutral scalar collapse that are plotted in 
Fig. (a) dynamical equation for r: rr^u~ — r^t, (b) for rj: rj^tt ~ (c) for a- ~ (d) for ip-. V’.tt ~ —2r^t'ip^tlr. 

(e) Inr = aln^ + 6, a = 0.5047 ± 0.0002, b = -0.286 ± 0.001. cr = aln^ + 6, a = -4.373 ± 0.002, b = -11.50 ± 0.01. (f) 
pj = a\n^ + b, a =-0.6087 ±0.0002, b =-1.131 ± 0.002. lnm = alnr±b, a =-18.307 ± 0.004, b =-31.39 ± 0.01. 
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(a) (b) 





Figure 17: (color online). Dynamics on the slice (x = 0.5, t = t) in null charge scattering. When the mass inflation happens, 
the held equations become null, in the sense that the temporal and spatial derivatives are almost equal. For example, in 
r,tt~r,xx and Near the center, (a) — rr,tt « —rr^xx~r‘ft~r%. (b) ~ rj^xx- (c) and (d): cr,tt « a^xx ~ ~ dTr 

(e) 'ip.tt ~ 'ij^.xx ~ ~ —‘2r^x'4>,xlr. (f) ratios between spatial and temporal derivatives for r, a, and ip. 
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•< -r_ 

r+ - J 

V A 

for r > r+ ^ 
for r < r+ 



0 0.5 1 1.5 2 


(e) 


Figure 18: Evolutions in null charge scattering on the slice (x = 0.5, t = t). (a)-(e): evolutions for r, a, tp, m, and |1 — K'^\. At 
the early stage of mass inflation, 1.2 < t < 2, r varies slowly; while at the late stage, t > 2, r varies faster and faster toward a 
small value close to zero. 
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(a) 






Figure 19: Evolutions at the late/fast stage of null charge scattering on the slice (x = 0.5,t = t). (a) Inr = aln^ + h, 
a = 0.4893 ± 0.0008, h = -0.956 ± 0.002. cr = aln^ + 6, a = -6.95 ± 0.05, h = 28.5 ±0.1. ip = aln^ ± 6, a = -0.717 ± 0.004, 
b = -0.24 ±0.01. (b) Inm = alnr ±6, a = -29.5 ± 0.2, b = 15.0 ± 0.4. 


inner horizon will decide which function each vari¬ 
able can take. On the other hand, we find that, as 
shown in Figs. [^a) and [^b) , the constraint equations 
@ and @ provide some useful information on the 
connections between r, a, and ip at the early stage of 
charge scattering: 


r,tcr,t 


4:TTr-ip\. 


(83) 


As shown in Fig. 18 the quantities r, a, ip, |1 — 
and m change dramatically at the beginning of charge 
scattering where Moreover, note that, near the 

central singularity, r, a, and ip have approximate analytic 
expressions in terms oi P, = Iq — t, where to is the time 
coordinate of the singularity curve. So it is natural to 
guess that, at the early stage of charge scattering, the 
above quantities may be also expressed by function of 
p = t — ts, where tg is a certain time value related to the 
early stage of charge scattering. We plot the evolutions 
of these quantities at the early stage of charge scattering 
in Fig. 1201 from which one can see that r, r^t, and cr may 
have the following approximate analytic expressions: 


r~r-, 




(84) 


(85) 


terms in the equation of motion for cr in charge scattering, 
at the very early stage (t « 1) of the collision between 
the scalar field and the inner horizon, compared to other 
terms, the contributions from the terms related to ip are 
tiny. Therefore, at this stage, the evolution of cr should 
not be much different from the corresponding one in the 
Reissner-Nordstrom geometry, 

As shown in Fig. 20’e), ip can be well fitted by power 
law functions of p, 


« cC ± IpQ. 


(87) 


We also plot (1 — K^) in Fig. 20 e), and find that ln(l — 
can be well fitted linearly with respect to 


Hl-K^)^fp + h. 


( 88 ) 


Currently we do not have derivations for this linear re¬ 
lation. The good thing is that, in the mass function, 
(1 — AT^) is a minor factor. Therefore, as verified in 
Fig. [20](f), the mass function can be reduced to 


r 

r_ 




. 2^7 2 


2 


(89) 


2b 


a^p^\ 


i 61 n^ ± CTo. 


( 86 ) 


In fact, a logarithmic expression for a is supported 
by its behavior near the inner horizon in the Reissner- 
Nordstrom black hole case. From Eqs. 0 and ( |2^ , one 
obtains that as r approaches the inner horizon r = r_, in 
the case of cr can be approximated by a logarithmic 
function of t. As shown in Fig. 0c), describing the 


where a, b, A, and erg are dehned in Eqs. (85) and (86). 
We list the fitting results below: 

(i) r t«a(f + by + d, a = (-2.38 ± 0.02) x b = 

-0.711 ± 0.001, c = 7.164 ± 0.007, d = (-2.395 ± 
0.004) X 10-4. 

(ii) afvaln{t + b) + c, a = 34.11 ± 0.02, 

b = -0.181 ± 0.001, c = 3.38 ± 0.03. 
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(a) 




t 



t 

(e) 



Figure 20: (color online). Numerical results along the slice (x = 0.5, t = t) in null charge scattering, (a) is for the {tx} constraint 
equation (181, and (b) is for the {tt} + {xx} one ( |19| l . (a) and (b) provide some useful information on the connections between 
some quantities at the early stage of charge scattering: ~ (c)-(f): evolutions at the early/slow stage of null charge 

scattering, (c) r,t = a{t + by + d, a = (-2.38±0.02) x lO'^ b = -0.786±0.001, c = 7.164±0.007, d = (-2.395±0.004) x 10"'‘. 
(d) (7 = aln(t + 6) + c, a = 34.11 ± 0.02, b = -0.256 ± 0.001, c = 3.38 ± 0.03. (e) i/i = a(t + b^ + d, a = (2.72 ± 0.06) x 10■^ 
b = -0.245 ± 0.004, c = 5.82 ± 0.01, d = (9.3 ± 0.1) x lO"'^. ln(l - K^) = at + b, a = -8.21 ± 0.02, b = 12.29 ± 0.04. (f) 
Inm = aln(t + 6)+ c, a = 77.52 ± 0.05, 6 =-0.201 ± 0.001, c =-16.15 ± 0.08. 
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(a) 



t 



r 

(c) 





r 

(d) 


Figure 21: (color online). Dynamics and solutions in critical charge scattering on the slice {x = 1.53,t = t). (a) and (b): 
dynamical equations for r and a. (c): evolutions of r and a. (d): evolutions of ip and m. 


(iii) + bY + d, a = (2.72 ± 0.06) x lO-^, b = 

-0.170 ± 0.004, c = 5.82 ± 0.01, d = (9.3 ± 0.1) x 
10 - 4 . 


(iv) \n{l-KY-at+b, a = -8.21±0.02, b = 11.67±0.03. 

(v) lnTO«oln(f + b) + c, a = 77.52 ± 0.05, b = —0.126 ± 
0.001, c= -16.16 ±0.08. 

It can be interesting to compare the dynamics at the 
early stage of null scattering with (1) that near the cent¬ 
ral singularity in spacelike scattering and (2) that at the 
late stage of null scattering. At the early stage of null 
scattering, r and Y change slowly: rrar-, r^t ~ with 
A > 1, and '0 « + ipo- Near the center, r and 0 

change fast: r « 0, x _ ^ ~ 1/2 < 

and ijjRiClnY In both cases, cr changes fast and has log¬ 
arithmic expressions with respect to 0 and respectively. 


E. Critical scattering 


In Secs. lynWvdi^ we studied spacelike and null 
scattering. In this subsection, we discuss critical scat¬ 
tering, which is on the edge between spacelike and null 
scattering. In fact, the dynamics of this type of mass 
inflation is similar to that at the late stage of null scat¬ 
tering. 


The dynamics of r, tr and the evolutions of r, cr, 0, 
and m in critical scattering on the slice {x = 1.53, t = t) 
are plotted in Fig. 21 In spacelike scattering, the terms 
in the field equations are grouped according to temporal 
and spatial derivatives. Take the field equation for a as 
an example. As shown in Fig. 16 c), at small-r regions, 
there are ~ 47r0^j and a^xx ~ ths early 

stage of the null scattering, the terms in the field equa¬ 
tions are grouped according to the types of quantities. As 
shown in Fig. c), at the early stage of the null scatter- 
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(a) 



(b) 




Figure 22: Evolutions in weak scalar charge scattering. The time interval between two consecutive slices is 30At = 0.15. (a) 
and (b): evolutions of r. In the current simulations, r does not approach zero. We expect this is also the case in later evolutions, 
(c) and (d): evolutions of cr and ip. 


ing, there are a^tt ~ o'^xx and In critical 

scattering, the dynamics has features from both spacelike 
and early stage of null scattering. As shown in Fig.[2^b), 
in this case, there are a,tt ~ cf,xx ~ ~ b’ur- 

ther details are skipped. 


VIII. WEAK SCALAR CHARGE SCATTERING 

In this section, we consider charge scattering with a 
weak scalar held. Parameter settings in this section are 
almost the same as those in the last section with the 
following exceptions: 

(i) Physical scalar held: 

ip{x,t)\t=o = aexp [-(a: - xqY/ h], 


a = 0.04, 6=1, and xq = 4. 

(ii) Grid. Spatial range: x e [—10 10]. Grid spacings: 
Ax = At = 0.005. 


As discussed in the above section, in some spacetime 
regions where the scalar held is strong, the inner horizon 
can contract to zero volume, and the central singularity 
becomes spacelike or null. However, this does not always 
necessarily happen. After all, it takes energy for the 
inner horizon to contract. When the scalar held carries 
less energy, the inner horizon may only contract to 
a nonzero value. This is conhrmed by our numerical 
results plotted in Fig. 22 which are in agreement with 
the numerical work in Ref. [29] and the mathematical 
proof in Ref. [S]. Since in this case the inner horizon is 
not totally destructed, one needs to reconsider whether 
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Figure 23: (color online). Dynamics and solutions in weak scalar charge scattering on the slice {x = 3,t = t). (a)-(d): dynamical 
equations for r, a, and ip. (e) and (f): evolutions of r and m. 





















































32 


the strong cosmic censorship conjecture is valid here. If 
not, more efforts need to be spent on understanding the 
nature of this conjecture and the internal structure of 
Reissner-Nordstrom and Kerr black holes directly. 

The dynamics for the quantities r, a, and tp are plotted 
in Figs. 23|a)-23 d). The numerical results show that at 
the later stage, the field equations for such quantities 
become null, in the sense that the temporal and spatial 
derivatives are almost equal, i.e., a^tt ~ o'^xx- Moreover, 
the derivatives have oscillations. As shown in Figs. 23'e) 
and[2^f), the mass function grows up continuously even 
as r approaches a constant value. Further details are 
skipped. 


IX. SUMMARY 

In this paper, we studied neutral scalar collapse, neut¬ 
ral scattering, and charge scattering numerically. Mass 
inflation was found to happen in the vicinity of the cent¬ 
ral singularity. Approximate analytic solutions were par¬ 
tially obtained. We summarize our work on computa¬ 
tional and physical issues separately below. 


A. Computational issues 


use this set of coordinates. For the {u, v) choice, in 
the field equations, many terms are mixed derivat¬ 
ives of u and v, e.g., r „„. However, for the {t,x) 
choice, spatial and temporal derivatives are usually 
separated, e.g., {r^tt - r^xx)- 

In the study of charge scattering, we set the initial 
conditions close to those in a Reissner-Nordstrom 
geometry. Consequently, removing the terms re¬ 
lated to the scalar field, we can test our code by 
comparing the numerical results to the analytic ones 
in a Reissner-Nordstrom geometry. Moreover, by 
comparing numerical results for charge scattering 
to the dynamics in a Reissner-Nordstrom geometry, 
we can obtain intuitions as to how the scalar field 
affects the geometry. 

(iv) Cauchy horizon: infinite or local regions? As im¬ 
plied by Eq. ( |11[ ), the exact inner horizon r = r- 
is at the regions where uv and (t^ — x^) are infin¬ 
ite. However, r still can be very close to the inner 
horizon even when uv and — x"^) take moder¬ 
ate values. Consequently, at regions where uv and 
— x^) take some moderate values, the scalar field 
and the inner horizon still can have strong interac¬ 
tions, resulting in mass inflation. 


(i) Numerical vs analytic approaches. Numerical and 
analytic approaches are both indispensable in grav¬ 
itational physics. The results from one of the two 
can be enlightening for the other one; and one ap¬ 
proach can be used to check the results from the 
other one. Typical examples include computations 
of gravitational waveforms in black hole binary sys¬ 
tems [5TH55] and studies of quantum evaporation of 
two-dimensional black holes [53]. The efficiency of 
combining the two approaches was displayed again 
in this paper. 

(ii) Problems with known solutions vs a new problem. 
In exploring charge scattering, we tried to closely 
compare problems with known solutions (dynamics 
in Schwarzschild and Reissner-Nordstrom geomet¬ 
ries and neutral scalar collapse) to a problem yet to 
be solved (charge scattering). 

(iii) dudv vs {—dt^ + dx'^) in double-null coordinates. In 
the studies of mass inflation, the dudv format of 
the Kruskal-like coordinates, ds^ = Ae~^'^dudv -\- 

is usually used. In this paper, we used the 
{—dt^ -\- dx'^) format instead, ds"^ = e“^'^(—-I- 
dx^) r^dVt^, with u = {t — x)j2 = const and 
V = {t x)/2 = const. Moreover, we set the initial 
conditions close to those in a Reissner-Nordstrom 
geometry. 

In the (t, a;) line element, one coordinate is time¬ 
like, and the rest are spacelike. We are used to this 
setup. It is more convenient and more intuitive to 


B. Physical issues 


(i) Interpretations of mass inflation in charge scatter¬ 
ing: analytic computations vs numerical simula¬ 
tions. According to the original papers on mass in¬ 
flation mma, the mechanism of mass inflation can 
be described below. Near the Cauchy horizon, the 
influx is infinitely blueshifted, while the outflux gen¬ 
erates the crucial separation between the Cauchy 
and inner horizons. Then the blueshift and redshift 
do not cancel, and the mass parameter increase dra¬ 
matically. In this paper, we explored mass inflation 
numerically. We examined the contributions of all 
the terms in the field equations and compared the 
dynamics in a Reissner-Nordstrom geometry and in 
charge scattering. In a Reissner-Nordstrom geo¬ 


metry, in the coordinates (101, r ^ asymptotes to 
zero near the inner horizon, and the mass parameter 
remains constant. In charge scattering, a scalar field 
impacts the inner horizon, and makes the metric 
quantity a in the line element grow faster than in 
the Reissner-Nordstrom geometry, causing a smaller 
repulsive force for r. Then r can cross the original 
inner horizon: mass inflation takes place. 

(ii) Contributions to mass inflation in charge scattering: 
e^'^ vs (r^ — r\). In a Reissner-Nordstrom black 
hole and charge scattering, the Misner-Sharp mass 
function can be written as 


r 

m = — 
2 




( 90 ) 
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Near the inner horizon of a Reissner-Nordstrom 
black hole, ^ 1. However, . This 

makes {r\ — r\) vanish. Then the mass function 
remains constant. Now we consider the early stage 
of null scattering that was discussed in Sec. |VIID 
We use the example plotted in Fig. At the 
slow stage, and (1 — vary roughly 

from 10® to 10^®, 10“^ to 10“^, and 1 to 10“^, re¬ 
spectively. So seems to play a more important 
role than r\ and (1 — r\/r\). However, compar¬ 
ing charge scattering and dynamics in the Reissner- 
Nordstrom metric, one can see that the mass func¬ 
tion increases dramatically not because is very 
large, but because r\ moves from extremely tiny 
values in a Reissner-Nordstrom geometry to a small 
number in charge scattering. Namely, r\ is more 
important than in the growth of the mass func¬ 
tion. 

(iii) Gravity vs repulsive (electric) force Nature of the 
central singularity (timelike, spacelike, or null). In 
charge scattering, the quantities in the field equa¬ 
tions can be separated into two sides: the gravitat¬ 
ing side (r, a, and ip) and the repulsive side (electric 
field). The dynamics in charge scattering mainly de¬ 
scribes how these quantities interact, including how 
the gravitating and repulsive sides compete. 

In the Reissner-Nordstrom geometry, at small-r re¬ 
gions, the electric field dominates gravity. As a res¬ 
ult, the central singularity is timelike. In charge 
scattering, when the scalar held is strong enough, 
the total gravity from the black hole and the scalar 
held dominates the electric held. The inner hori¬ 
zon contracts to zero volume rapidly, and the cent¬ 
ral singularity becomes spacelike. When the scalar 
held is intermediate, the inner horizon contracts to 
a small or zero value. The equations of motion be¬ 
come null, in the sense that in the equations, the 
temporal and spatial derivatives are almost equal, 
e.g., cr tt « cr^xx- In the case of the inner horizon 
contracting to zero, the central singularity becomes 
null. 

(iv) Compare neutral collapse and different types of 
charge scattering. In the late stages of strong and 
intermediate scalar charge scattering, the dynamics 
and hence the solutions in the vicinity of the cent¬ 
ral singularity are similar to those in neutral scalar 
collapse. 

Regarding the early stage of intermediate scalar 
charge scattering, we noticed that, in the mass func¬ 
tion, is the dominant factor, and hence the mass 
function can be well fitted by a logarithmic function 
of the coordinate time t. 

(v) The inner horizon in a Reissner-Nordstrom black 
hole vs the central singularity in a Schwarzschild 
black hole. These two share some similarities as 
below. 


For Reissner-Nordstrom and Schwarzschild black 
holes, throughout the whole spacetime, the Misner- 
Sharp mass function is constant. When a scalar field 
impacts the inner horizon of a Reissner-Nordstrom 
black hole, the scalar field can modify the geometry 
in the vicinity of the inner horizon significantly, es¬ 
pecially on r t. The inner horizon contracts and 
mass inflation takes place. In neutral scalar col¬ 
lapse toward a Schwarzschild black hole formation, 
the scalar field can also modify the geometry in the 
vicinity of the central singularity dramatically, es¬ 
pecially on the metric component cr |20j . Then mass 
inflation also happens. 

The BKL conjecture is an important result on the 
dynamics in the vicinity of a spacelike singular¬ 
ity [T 2 HI 5 ]. The first statement of this conjecture 
is that as the singularity is approached, the dy¬ 
namical terms dominate the spatial terms in the 
field equations. In other words, the way gravity 
changes over time is more important than the vari¬ 
ation of the gravitational held from one location 
to the next m- We would like to say that, to 
a large extent, later evolutions in a strong gravit¬ 
ational held largely erase away the initial informa¬ 
tion on the connections between neighboring points. 
As discussed in Ref. [5D] and also in this paper, in 
double-null coordinates, using the above arguments, 
one can interpret the following behaviors displayed 
in numerical simulations: near the central singu¬ 
larity of a Schwarzschild black hole and the inner 
horizon of a Reissner-Nordstrom black hole, there 
is 


f’.x 


r,x 

V'.t 




(91) 


In this paper, it was found that Eq. (911 can inter¬ 
pret how mass inhation takes place. 


The second and third statements of the BKL con¬ 
jecture are that i) the metric terms dominate the 
matter held terms, while the matter held may not 
be negligible if it is a scalar held; ii) the dynam¬ 
ics of the metric components and the matter helds 
is described by the Kasner solution. These two 
statements were conhrmed in simulations of neutral 
scalar collapse in f{R) gravity in Ref. [50] and in 
general relativity in this paper. The second state¬ 
ment was also verihed in charge scattering in this 
paper. However, the third statement on Kasner 
solution may not apply to charge scattering. 


Near the central singularity in a Schwarzschild 
black hole, the asymptotic solutions are r«A^ 2 , 
CT«Hln^ -I- (To, and ippeChi(, where £,=to — t, and 
to is the time coordinate on the singularity curve 
r{x, t) = 0. Near the inner horizon of a Reissner- 
Nordstrom black hole, the asymptotic solutions are 
-be with 6 > 1, cr«dln^-|-/, and ip^hf^ +'0Oj 
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where ^=t — ts, and ts is the coordinate time t at 
the beginning of charge scattering. 

(vi) Strong cosmic censorship conjecture. On the 
Cauchy horizon of a Reissner-Nordstrom black hole, 
the predicability is violated. Based on first-order 
calculation, the inner (Cauchy) horizon appears un¬ 
stable under perturbations [7]. Then it was con¬ 
jectured that for generic asymptotically flat initial 
data, the maximal Cauchy development is future in¬ 
extendible. In this paper, it was found that for weak 
scalar field perturbation, the inner (Cauchy) hori¬ 
zon only contracts to a nonzero value. This is also 
in agreement with the numerical work in Ref. |29j 
and the mathematical proof in Ref. [9]. Therefore, 
one needs to reconsider the validity of the strong 
cosmic censorship. 


For reference concern, we list the derivations here via 
analog to the Schwarzschild black hole case in the form 
of Ref. [12] . 

We start from the Reissner-Nordstrom metric in the 
conventional form, 


ds^ = -Adf + A-^dr^ + r^dVt^, 


where 


A _ 1 ‘2.m _{r- r+)(r - r_) 

^ ^ ' 2 2 ’ 


(Al) 


(A2) 


= m ± \/ m? — q^. (A3) 

Considering the radial null curves along which = 0, 
we have 


(vii) Inside vs outside blaek holes: loeal vs global. 
Throughout the whole spacetime of stationary 
Schwarzschild and Reissner-Nordstrom black holes, 
the Misner-Sharp mass function is equal to the black 
hole mass. For a gravitational collapsing system, at 
asymptotic flat regions, the mass function describes 
the total mass of the dynamical system. However, 
in this system, near the central singularity of a 
Schwarzschild black hole or near the inner horizon 
of a Reissner-Nordstrom black hole, the dynamics 
is local. Then the mass function does not provide 
global information on the mass of the collapsing sys¬ 
tem. 

In summary, in this paper, we studied neutral scalar 
collapse, neutral scattering, and charge scattering. Re¬ 
garding charge scattering, for convenience and intuitive¬ 
ness, Kruskal-like coordinates were used, and initial con¬ 
ditions were set up to be close to those in a Reissner- 
Nordstrom geometry. Mass inflation was also found to 
happen near the central singularity in neutral scalar col¬ 
lapse. Approximate analytic solutions for mass inflation 
were partially obtained. Connections between Schwarz¬ 
schild black holes, Reissner-Nordstrom black holes, neut¬ 
ral collapse, neutral scattering, and charge scattering 
were explored. 


dr \ r 


(A4) 


The solution to the above equation can be expressed as 
t = ±r* -I- const, (A5) 

where 

(A6) 

(A7) 


1 , 

r 

1 , 

r 

r + —— In 

1- 

+ — In 

1- 

2fc_|_ 

r+ 

2k- 

r_ 


k+ = 


r± - r:^ 

24 ' 


Define 


u = t — r, v = t + r 


and then the line element (Al I becomes 


= A{—dr + dr*'^) -|- r^dil^ 

= — ^A{dudv + dvdu) + r'^dlt^. 


(AS) 


(A9) 


1. Patch I: r_ < r < +00 
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Appendix A: Reissner-Nordstrom metric in 
Kruskal-like coordinates 

The Reissner-Nordstrom metric in Kruskal-like co¬ 
ordinates was obtained in Ref. |3T| (also see Ref. [12]) ■ 


(i) r+ < r < - 1-00 of Patc h I 
As implied in Eq. (A6), the inner and outer horizons 
in the metric (A9) are pushed to infinity. To compact 
them into finite space, we define 

I/+=-6 -'^'+“, U+=e'=+^. for r>r+. (AlO) 


Rewrite Eq. (A61 as 


1 

1 


In 


r 

/ r \ 

1 - 

-1 

r+ 

J\ 


1-h 


\k-\ 


In 


r 

1 - 

r_ 


(All) 
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Using Eqs. (A2), (A8), (AlOl, and (All), the line ele¬ 


ment (A9) can be converted into the following format 

k. 


ds^ = - 


r+r_ 2 fe 


\r_ 


-^ 


1 + 




(A12) 

Note that we get used to coordinate systems in which 
one coordinate is timelike and the rest are spacelike. So 
defining 


T = 


+ U+ . 


= 6*^+’’ sinh(fc+t), 


(A13) 


In a similar routine as in Patch 1, the metric for the 
spacetime of r_ < r < r+ can be obtained as 


ds^ = ("l - {-dT^+dX^)+r^dn\ 

k_r^ \ r+ / 


with 


r2_x2=e-2|'=-l’-f--lVl--^ 

r_ J \ r+ 


(A23) 


(A24) 


U+ - t/+ t. . 

X = --- = e^+’' cosh(fc+t), 


(A14) 


one can rewrite the metric (A12) as 

fc I 




k\r'^ 


r_ 


i+- 


*-i 


Moreover, using Eqs. (A6), (A13), and (A14), r and t 
can be expressed as 


{-dT^+dX^)+r^dn\ 

(A15) 


T^-X^ = e2'=+'' (1--1 

r+J Vr_ 


|fc-i 


(A16) 


— = coth(A:_<). (A25) 

JC 

(ii) 0 < r < r_ of Patch II 
Define 

[/- = -e-'=-“, U-=e'=-U forO<r<r_. (A26) 

Then we have 


T 

— = tanh(fc_t). 


(A27) 


The line element is the same as Eq. (A23 l, and the ex¬ 
pression for r is the same as Eq. (A241. 


T 

— = tanh(A:+t). 

y(. 

(ii) T- < r < r+ of Patch I 
Define 


(A17) 


U+ = -e-'=+“, U+ = -e'^+U for r_ < r < r+. 

(A18) 

Similar to the case of < r < -|-oo, one can obtain 


T = 


+ U+ 


= cosh(fc+t), (A19) 


U+ - 1/+ i, . 

X = --- = sinh(A:+f), 


T 


— = coth(fc_|_t). 
X 


(A20) 


(A21) 


The line element is the same as Eq. (A15l, and the ex¬ 


pression for r is the same as Eq. (A16l 


2. Patch II: 0 < r < r+ 


Appendix B: Einstein tensor and Energy-momentum 
tensor of a massive scalar field 


In this appendix, we give specific expressions of the 
Einstein tensor and the energy-momentum tensor of a 
massive scalar field. In double-null coordinates expressed 
by Eqs. ([^ and (§, some components of the Einstein 
tensor can be expressed as follows: 


Gl = 


2e 


la r 


* ^2 


r{r,t<y,t + r^x<y,x) +rrx 




G X _ 

^ — 


2 e' 


2a r 


- + r^x(J,x) - rru 




(Bl) 

(B2) 

(B3) 

(B4) 


Gq — G^ — + T^xx — + cr^xx )]) (B5) 


(i) r_ < r < r+ of Patch II 
Define 


U- = -e 


— k—u 


V = —for r_ < r < r+. 


(A22) 


Guu — i'^,uu 
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Gyy - (j’,VV “ 1 “ 2(7 
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(B6) 
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For a massive scalar field with energy-momentum 
tensor 


there are 




Qfj-iy 


T/ = -e^'" 




rpO _ rp<p _ 2(7 

J- f) — J- ^ ^ 


7,0^1+ ^%)+e-^'^VW 






Tuu — 


Tvv = 


(B8) 


, (B9) 


(BIO) 


(Bll) 


(B12) 


(B13) 




The equations obtained in this appendix can be used to 
derive the field equations as discussed in Sec. |III A| 


Appendix C: Schwarzschild geometry in Kruskal 
coordinates 


Y can be a negative or a complex number. 

The Hrst- and second-order derivatives of W are 


dW 

dz 

d?W _ 
dz"^ 


W 

z(l + W)’ 

W^(2 + W) 

'z^(l + W)3^ 


for zj^iO,- 


1 


for z ^ < 0, — 


1 


(C6) 


(C7) 


Consequently, with Eqs. (|C4|, (C6), and (C7|, one ob¬ 


tains the first- and second-order derivatives of r with re¬ 
spect to x: 

1 dr dW 2x 
2m dx dz e ’ 


(C8) 


1 d'^r 
2m dx^ 


~d^ 


2a;y dW 
e ) dz e 


(C9) 


Near the singularity curve, z[= (a;^ — t^)/e] approaches 
— 1/e, and W asymptotes to —1. Consequently, the 
second-order derivative of r with respect to x can be ap¬ 
proximated as follows: 


(B14) _L . ~ _ 

dx^ 


1 

2m 


4x^ 


d^W (2x 


(1 -I- WY dz'^ \ e 


(CIO) 


Similarly, one obtains the first- and second-order deriv¬ 
atives of r with respect to t near the singularity curve: 


1 

2m 


dr 

dt 


dW 

dz 


2t 


(Cll) 


In this appendix, we derive the analytic expressions for 
the spatial and temporal derivatives near the singularity 
curve for a Schwarzschild black hole in Kruskal coordin¬ 


ates, and verify the mass formula (46). The Schwarz 


schild metric in Kruskal coordinates is 

d^2 = 32m e-^(^-df dx^) + r'^dVi^, (Cl) 

r 


with 


2_a;2= 

V 2m/ 


P ■= - P 2 

32m^ 


The solution to Eq. (C2) is 


2m 


= 1 + W{z), 


where 


x^-P 


z = 


and W is the Lambert W function defined by 
Y = W{Y)e^(^\ 


(C2) 


(C3) 


(C4) 


(C5) 


1 

2m 


dr 

dt^ 


4C 


d^W 


(1 + VF)3 dz^ ye J ' 

Therefore, with Eqs. (C8) and ( |C10 )-(C12), the ratios 
between the spatial and temporal derivatives can be ex¬ 
pressed by the slope of the singularity curve, K, 


r x X 

-^ = -= -K, 

r,t t 


• ,XX 

r,tt 


(f) 


- =K-^ 


(C13) 


(C14) 


As discussed in Sec. EYl in spherical scalar collapse in 
double-null coordinates, in the vicinity of the singularity 
curve of the formed black hole, for the metric compon¬ 
ents and scalar field, the ratios between the spatial and 
temporal derivatives are also defined by K. 


Using Eqs. ( |C2[ ), ( |C3[ ), ( |C4[ ), ( |C8[ ), and ( |C11[ ), one ob¬ 
tains 


rl-r^ =64m4 


(a;2 — t'^)r^ 


32m3 




(CIS) 
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This implies that the definition of mass function (461 for 
a Schwarzschild black hole is valid throughout the space- 
time of the black hole, inside and outside the horizon, 
and also in the vicinity of the central singularity, 


In other words, the Misner-Sharp mass function is equal 
to the black hole mass everywhere. 
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